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The blackbody radiation left over from the Big Bang has been transformed by 
the expansion of the Universe into the nearly isotropic 2.73 K Cosmic Microwave 
Background. Tiny inhomogeneities in the early Universe left their imprint on the 
microwave background in the form of small anisotropies in its temperature. These 
anisotropics contain information about basic cosmological parameters, particularly 
the total energy density and curvature of the universe. Here we report the first images 
of resolved structure in the microwave background anisotropies over a significant part 
of the sky. Maps at four frequencies clearly distinguish the microwave background 
from foreground emission. We compute the angular power spectrum of the microwave 
background, and find a peak at Legendre multipole ipeak = (197 it 6), with an amplitude 
^^200 = (69 ± 8)fj,K. This is consistent with that expected for cold dark matter models 
in a flat (euclidean) Universe, as favoured by standard inflationary scenarios. 

Photons in the early Universe were tightly coupled to ionized matter through Thomson 
scattering. This coupling ceased about 300,000 years after the Big Bang, when the Universe 
cooled sufficiently to form neutral hydrogen. Since then, the primordial photons have travelled 
freely through the universe, redshifting to microwave frequencies as the universe expanded. We 
observe those photons today as the cosmic microwave background (CMB). An image of the early 
Universe remains imprinted in the temperature anisotropy of the CMB. Anisotropies on angular 
scales larger than ~ 2° are dominated by the gravitational redshift the photons undergo as they 
leave the density fluctuations present at decoupling^^'^) . Anisotropies on smaller angular scales are 
enhanced by oscillations of the photon-baryon fluid before decoupling(^) . These oscillations are 
driven by the primordial density fluctuations, and their nature depends on the matter content of 
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the universe. 

In a spherical harmonic expansion of the CMB temperature field, the angular power spectrum 
specifies the contributions to the fluctuations on the sky coming from different multipoles, each 
corresponding to the angular scale = 7r/£. Density fluctuations over spatial scales comparable to 
the acoustic horizon at decoupling produce a peak in the angular power spectrum of the CMB, 
occurring at multipole ipeak- The exact value of £peak depends on both the linear size of the 
acoustic horizon and on the angular diameter distance from the observer to decoupling. Both 
these quantities are sensitive to a number of cosmological parameters (see for example ref.4), but 
ipeak primarily depends on the total density of the Universe, fio • In rnodels with a density 
near 1, ipeak ~ 200/0^/^. A precise measurement of peak can efHciently constrain the density and 
thus the curvature of the Universe. Observations of CMB anisotropies require extremely sensitive 
and stable instruments. The DMR^ instrument on the COBE satellite mapped the sky with an 
angular resolution of ~ 7°, yielding measurements of the angular power spectrum at multipoles 
i < 20. Since then, experiments with finer angular resolution^" have detected CMB fluctuations 
on smaller scales and have produced evidence for the presence of a peak in the angular power 
spectrum at peak ipeak ~ 200. 

Here we present high resolution, high signal-to-noise maps of the CMB over a significant fraction 
of the sky, and derive the angular power spectrum of the CMB from i= 50 to 600. This power 
spectrum is dominated by a peak at multipole ipeak = (197 =•= 6) (la error). The existence of this 
peak strongly supports inflationary models for the early universe, and is consistent with a flat, 
Euclidean Universe. 

The Instrumrent 

The BOOMERanG (Balloon Observations Of Millimetric Extragalactic Radiation and 
Geomagnetics) experiment is a microwave telescope that is carried to an altitude of ~ 38 km by 
a balloon. BOOMERanG combines the high sensitivity and broad frequency coverage pioneered 
by an earlier generation of balloon-borne experiments with the long (~ 10 days) integration time 
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available in a long-duration balloon flight over Antarctica. The data described here were obtained 
with a focal plane array of 16 bolometric detectors cooled to 0.3 K. Single-mode feedhorns provide 
two 18' full-width at half-maximum (FWHM) beams at 90 GHz and two 10' (FWHM) beams at 
150 GHz. Four multi-band photometers each provide a 10.5, 14 and 13 FWHM beam at 150, 240 
and 400 GHz respectively. The average in-flight sensitivity to CMB anisotropies was 140, 170, 210 
and 2700 /j^K ■ s^/^ at 90, 150, 240 and 400 GHz respectively. The entire optical system is heavily 
baffled against terrestrial radiation. Large sun-shields improve rejection of radiation from > 60° in 
azimuth from the telescope borcsight. The rejection has been measured to be greater than 80 dB 
at all angles occupied by the Sun during the CMB observations. Further details on the instrument 
can be found in refs 17-21. 

Observation 

BOOMERanG was launched from McMurdo Station (Antarctica) on 29 December 1998, at 
3:30 GMT. Observations began 3 hours later, and continued uninterrupted during the 259-hour 
flight. The pay load approximately followed the 79° S parallel at an altitude that varied daily 
between 37 and 38.5 km, returning within 50 km of the launch site. 

We concentrated our observations on a target region, centred at roughly right ascension (RA) 
5h, declination (dec.) —45°, that is uniquely free of contamination by thermal emission from 
interstellar dust*^^^^ and that is approximately opposite the Sun during the austral summer. We 
mapped this region by repeatedly scanning the telescope through 60° at fixed elevation and at 
constant speed. Two scan speeds (l°-s~-'^ and 2°-s~-'^ in azimuth) were used to facilitate tests 
for systematic effects. As the telescope scanned, degree-scale variations in the CMB generated 
sub-audio frequency signals in the output of the detector^^^^ . The stability of the detetector system 
was sufficient to allow sensitive measurements on angular scales up to tens of degrees on the sky. 
The scan speed was sufficiently rapid with respect to sky rotation that identical structures were 
observed by detectors in the same row in each scan. Detectors in different rows observed the same 
structures delayed in time by a few minutes. 
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At intervals of several hours, the telescope elevation was interchanged between 40°, 45° and 50° 
in order to increase the sky coverage and to provide further systematic tests. Sky rotation caused 
the scan centre to move and the scan direction to rotate on the celestial sphere. A map from a 
single day at a single elevation covered roughly 22° in declination and contained scans rotated by 
±11° on the sky, providing a cross-linked scan pattern. Over most of the region mapped, each 
sky pixel was observed many times on different days, both at l°-s~-'^ and 2°-s~^ scan speed, with 
different topography, solar elongation and atmospheric conditions, allowing strong tests for any 
contaminating signal not fixed on the celestial sphere. 

The pointing of the telescope has been reconstructed with an accuracy of 2' r.m. s.using data 
from a Sun sensor and rate gyros. This precision has been confirmed by analyzing the observed 
positions of bright compact HIT regions in the Galactic plane (RCWSS^^ , RCW57, IRAS08576 
and IRAS1022) and of radio-bright point sources visible in the target region (the QSO 0483-436, 
the BL-Lac object 0521-365 and the blazar 0537-441). 

Calibrations 

The beam pattern for each detector was mapped before flight using a thermal source. The 
main lobe at 90, 150 and 400 GHz is accurately modelled by a Gaussian function. The 240 GHz 
beams are well modelled by a combination of two Gaussians. The beams have small shoulders 
(less than 1% of the total solid angle), due to aberrations in the optical system. The beam-widths 
were confirmed in flight via observations of compact sources. By fitting radial profiles to these 
sources we determine the effective angular resolution, which includes the physical beamwidth 
and the effects of the 2' r.m. s.pointing jitter. The effective FWHM angular resoluti on of the 
150 GHz data that we use here to calculate the CMB power spectrum is (lOitl)', where the error 
is dominated by uncertainty in the pointing jitter. 

We calibrated the 90, 150 and 240 GHz channels from their measured response to the CMB dipole. 
The dipole anisotropy has been accurately (0.7%) measured by COBE-DMR^^, fills the beam and 
has the same spectrum as the CMB anisotropics at smaller angular scales, making it the ideal 
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calibrator for CMB experiments. The dipole signal is typically ~ 3 mK peak-topeak in each 60° 
scan, much larger than the detector noise, and appears in the output of the detectors at f=0.008Hz 
and f=0.016Hz in the l°-s-^ and 2°-s-^ scan speeds, respectively. The accuracy of the calibration 
is dominated by two systematic effects: uncertainties in the low-frequency transfer function of the 
electronics, and low-frequency, scan-synchronous signals. Each of these is significantly different at 
the two scan speeds. We found that the dipole-fitted amplitudes derived from separate analysis 
of the l°-s~^ and 2°-s~^ data agree to within ±10% for every channel, and thus we assign a 10% 
uncertainty in the absolute calibration. 

From detector signals to CMB maps 

The time-ordered data comprises 5.4 • 10^ 16-bit samples for each channel. These data 
are flagged for cosmic-ray events, elevation changes, focal-plane temperature instabilities, and 
electromagnetic interference events. In general, about 5% of the data for each channel are flagged 
and not used in the subsequent analysis. The gaps resulting from this editing are filled with a 
constrained realization of noise in order to minimize their effect in the subsequent filtering of the 
data. The data are deconvolved by the bolometer and electronics transfer functions to recover 
uniform gain at all frequencies. 

The noise power spectrum of the data and the maximum-likelihood maps^^~^^ are calculated 
using an iterative technique^^ that separates the sky signal from the noise in the time-ordered 
data. In this process, the statistical weights of frequencies corresponding to angular scales larger 
than 10° on the sky are set to zero to filter out the largest-scale modes of the map. The maps are 
pixelized according to the HEALPix pixelization scheme^'^. 

Figure 1 shows the maps obtained in this way at each of the four frequencies. The 400 GHz 
map is dominated by emission from interstellar dust that is well correlated with that observed 
by the IRAS and COBE/DIRBE satellites. The 90, 150 and 240 GHz maps are dominated by 
degree-scale structures that are resolved with high signal-to-noise ratio. A qualitative but powerful 
test of the hypothesis that these structures are CMB anisotropy is provided by subtracting one 
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map from another. The structures evident in all three maps disappear in both the 90—150 GHz 
difference and in the 240—150 GHz difference, as expected for emission that has the same spectrum 
as the CMB dipole anisotropy used to calibrate the maps. 

To quantify this conclusion, we performed a "colour index" analysis of our data. We selected 
the ~ 18000 14' pixel at galactic latitude b < —15°, and made scatter plots of 90 GHz versus 
150 GHz and 240 GHz versus 150 GHz. A linear fit to these scatter plots gives slopes of 1.00 it 0.15 
and 1.10 ±0.16, respectively (including our present 10% calibration error), consistent with a CMB 
spectrum. For comparison, free-free emission with spectral index —2.35 would produce slopes of 
2.3 and 0.85, and is therefore rejected with > 99% confidence; emission from interstellar dust with 
temperature = 15 K and spectral index of emissivity a = 1 would produce slopes of 0.40 and 
2.9. For any combination of > 7 K and 1 < a < 2, the dust hypothesis is rejected with > 99% 
confidence. We conclude that the dominant source of structure that we detect at 90, 150 and 
240 GHz is CMB anisotropy. 

We further argue that the 150 GHz map at 6 < —15° is free of significant contamination by 
any known astrophysical foreground. Galactic synchrotron and free-free emission is negligible at 
this frequency^^. Contamination from extra-galactic point sources is also small^^; extrapolation 
of fluxes from the PMN survey^^ limits the contribution by point sources (including the three 
above-mentioned radio-bright sources) to the angular power spectrum derived below to < 0.7% 
at i = 200 and < 20% at i = 600. The astrophysical foreground that is expected to dominate at 
150 GHz is thermal emission from interstellar dust. We placed a quantitative limit on this source 
of contamination as follows. We assumed that dust properties are similar at high (6 < —20°) 
and moderate (—20° < 6 < — 5°) Galactic latitudes. We selected the pixels at moderate Galactic 
latitudes and correlated the structure observed in each of our four bands with the IRAS/DIRBE 
map, which is dominated by dust in cirrus clouds. The best-fit slope of each of the scatter plots 
measures the ratios of the dust signal in the BOOMERanG channels to the dust signal in the 
IRAS/DIRBE map. We found that the 400 GHz map is very well correlated to the IRAS/DIRBE 
map, and that dust at 6 < —20°) can account for at most 10% of the signal variance at 240 GHz, 
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3% at 150 GHz and 0.5% at 90 GHz. 
Angular Power Spectra 

We compared the angular power spectrum of structures evident in Fig.l with theoretical 
predictions. In doing so, we separated and removed the power due to statistical noise and 
systematic artifacts from the power due to CMB anisotropies in the maps. The maximum- 
likelihood angular power spectrum of the maps was computed using the MADCAP^^ software 
package, whose algorithms fully take into account receiver noise and filtering. 

Pull analysis of our entire data set is under way. Because of the computational intensity of 
this process, we report here the results of a complete analysis of a limited portion of the data 
chosen as follows. We analysed the most sensitive of the 150 GHz detectors. We restricted the 
sky coverage to an area with RA > 70°, b < —20° and —55° < Dec < —35°, and we used only 
the ~ 50% of the data from this detector that was obtained at a scan speed of l°-s~^. We used 
a relatively coarse pixelization of 8,000 14-arcmin pixels as a compromise between computation 
speed and coverage of high multipoles. Pinally, we limited our analysis to £ < 600 for which the 
effects of pixel shape and size and our present uncertainty in the beam size (1') are small and can 
be accurately modelled. 

The angular power spectrum determined in this way is shown in Pig. 2 and reported in 
Table 1. The power spectrum is dominated by a peak at ipeak ~ 200, as predicted by inflationary 
cold dark matter models. These models additionally predict the presence of secondary peaks. 
The data at high £ limit the amplitude, but do not exclude the presence, of a secondary peak. 
The errors in the angular power spectrum are dominated at low multipoles {£ < 350) by the 
cosmic/sampling variance, and at higher multipoles by detector noise. 

The CMB angular power spectrum shown in Pig.2 was derived from 4.1 days of observation. 
As a test of the stability of the result, we made independent maps from the first and second halves 
of these data. The payload travels several hundred kilometers, and the Sun moves 2° on the sky. 
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between these maps. Comparing them provides a stringent test for contamination from sidelobe 
pickup and thermal effects. The angular power spectrum calculated for the difference map is 
shown in Fig. 2. The reduced of this power spectrum with respect to zero signal is 1.11 (12 
degrees of freedom) , indicating that the difference map is consistent with zero contamination. 

A peak at £ ^ 200 implies a flat Universe 

The location of the first peak in the angular power spectrum of the CMB is well measured by 
by this data setthe BOOMERanG data set. Prom a parabolic fit to the data at 1=50 to 300 in the 
angular power spectrum, we find ipeak = (197 it 6) (Icr error). The parabolic fit does not bias the 
determination of the peak molutpolemultipole: applying this method to Monte Carlo realizations 
of theoretical power spectra we recover the correct peak location for a variety of cosmological 
models. Finally, Notice the peak location is independent of the details of the data calibration, 
which obviously affect only the height of the peak and not its location. The height of the peak is 
^^200 = (69 ± 4) ± 7 mK [l-a statistical and calibration errors, respectively). 

The data are inconsistent with current models based on topological defects (see, for example, 
ref. 35) but are consistent with a subset of cold dark matter models. We generated a database 
of cold dark matter models^^''^'', varying six cosmological parameters (the range of variation 
is given in parentheses): the non relativistic matter density, ^lm (0.05-2); the cosmological 
constant, f^A (0-1); the Hubble constant, h (0.5-0.8); the baryon density, h^O;, (0.013-0.025), the 
primordial scalar spectral index, ns (0.8-1.3); and the overall normalization A (freeparameter) of 
the primordial density fluctuation power spectrum. We compared these models with the power 
spectrum we report here to place constraints on allowed regions in this 6-parameter space. In 
Figure 3 we mark with black dots the region of the flm — plane where some combination of the 
remaining four parameters within the ranges defined by our model space gives a power spectrum 
consistent with our 95% confidence interval for ipeak- This region is quite narrow and elongated 
along the "flat Universe" line Um + ^^A = 1- The width of this region is determined by degeneracy 
in the models, which produce closely similar spectra for different values of the parameters'^*^. We 
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further evaluated the hkeUhood of the models given the BOOMERanG measurement and the 
same priors (constraints on the values of the cosmological parameters) as in ref. 16. Marginalizing 
over all the other parameters, we found the following 95% confidence interval for 0,q = + 
0.88 < < 1.12. This provides evidence for a euclidean geometry of the Universe. Our data 
clearly show the presence of power beyond the peak at ^ = 197, corresponding to smaller-scale 
structures. The consequences of this fact will be fully analyzed elsewhere. 
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TABLE 1: Angular power specrtum of CMB anisotropy 



^— range 


150 GHz 


150 GHz 




([1^* half] + [2""' half])/2 


([1** half]- [2"'^ half])/2 


[26-75] 


1140±280 


63±32 


[76-125] 


3110±490 


16±20 


[126-175] 


4160il=540 


17±28 


[176-225] 


4700±540 


59±44 


[226-275] 


4300±460 


68±59 


[276-325] 


2640±310 


130±82 


[326-375] 


1550±220 


-7±92 


[376-425] 


1310±220 


-60±120 


[426-475] 


1360±250 


0±160 


[476-525] 


1440±290 


220±230 


[526-575] 


1750±370 


130±300 


[576-625] 


1540±430 


-430±360 



Shown are measurements of the angular power spectrum of the cosmic microwave background at 
150 GHz, and test for systematic effects. The values listed arc for I\Tj = £[i + 1)^/271, in mK^. 
Here q =< aj^ >, and < > are the coefficients of the spherical harmonic decomposition of 
the CMB temperature field: AT{9,(f)) = ^ aemXemi^ , <!>) ■ The stated la errors include statistical 
and cosmic/sample variance, and do not include a 10% calibration uncertainty. 
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Figure 1: BOOMERanG sky maps (equatorial coordinates). The sky maps at 90, 150, and 
240 GHz (left panels) are shown with a common color scale, using a thermodynamic temperature 
scale chosen such that CMB anisotropics will have the same amplitude in the three maps. Only the 
colour scale of the 400 GHz map (bottom right) is 14 times larger than the others: this has been 
done to facilitate comparison of emission from interstellar dust (ISD), which dominates this map, 
with ISD emission present in the lower-frequency maps. The maps at 90 and 400 GHz are each from 
a single detector, while maps at 150 and 240 GHz have each been obtained by co-adding data from 
three detectors. For purposes of presentation, the maps have been smoothed with gaussian filters 
to obtain FWHM effective resolution of 22.5' (small circle in the bottom right side of each panel). 
Structures along the scan direction larger than 10° are not present in the maps. Several features are 
immediately evident. Most strikingly, the maps at 90, 150, and 240 GHz are dominated by degree- 
scale structures that fill the map, have well-correlated morphology and are identical in amplitude 
in all three maps. These structures arc not visible at 400 GHz. The 400 GHz map is dominated 
by diffuse emission which is correlated with the ISD emission mapped by IRAS/DIRBE^^. This 
emission is strongly concentrated towards the right-hand edge of the maps, near the plane of the 
Galaxy. The same structures are evident in the 90, 150 and 240 GHz maps at galactic latitude 
b > —15°, albeit with an amplitude that decreases steeply with decreasing frequency. 
The large-scale gradient evident especially near the right edge of the 240 GHz map is a result 
of high-pass filtering the very large signals near the Galactic plane (not shown). This effect is 
negligible in the rest of the map. The two top right panels show maps constructed by differencing 
the 150 and 90 GHz maps and the 240 and 150 GHz maps. The difference maps contain none of 
the structures that dominate the maps at 90, 150 and 240 GHz, indicating that these structures do 
indeed have the ratios of brightness that are unique to the CMB. The morphology of the residual 
structures in the 240 — 150 GHz map is well-correlated with the 400 GHz map, as is expected if 
the residuals are due to the ISD emission. Three compact sources of emission are visible in the 
lower-frequency maps, as indicated by the circles. These are known radio-bright quasars from the 
BEST pointing catalogue at 230 GHz. The boxed area has been used for computing the angular 
power spectrum shown in Fig.2. 
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Figure 2: Angular power spectrum measured by BOOMERanG at 150 GHz. Each point is the 
power averaged over A£ = 50 and has neghgible correlations with the adjacent points. The error 
bars indicate the uncertainty due to noise and cosmic/sampling variance. Tlic errors are dominated 
by cosmic/sampling variance at £ < 350; they grow at large I due to the signal attenuation caused 
by the combined effects^^ of the 10' beam and the 14' pixelization (0.87 at £=200 and 0.33 at 
£=600). The current ±10% uncertainty in the calibration corresponds to an overall re-scaling of 
the y-axis by ± + 20%, and is not shown. The current 1' uncertainty in the angular resolution of 
the measurement creates an additional uncertainty, indicated by the distance between the ends of 
the red error bars and the blue horizontal lines, that is completely correlated and is largest (11%) 
at i = 600. The green points show the power spectrum of a difference map obtained dividing the 
data in two parts corresponding to the first and second halves of the timestream. We make two 
maps (A and B) from these halves, and the green points show the power spectrum computed from 
the difference map, (A-B)/2. Signals originating from the sky should disappear in this map, so this 
is a test for contamination in the data (see text). The solid curve has parameters {Cli,, Jl^, Qa, 
ns, h) = (0.05, 0.31, 0.75, 0.95, 0.70) . It is the best fit model for the BOOMERanG test flight 
data^^'^^, and is shown for comparison only. The model that best fits the new data reported here 
will be presented elsewhere. 



Figure 3: Observational constraints on Q^n and ^a. All the cosmological models (from our 
data base) consistent with the position of the peak in the angular power spectrum measured by 
BOOMERanG (95% confidence intervals) define an "allowed" region in the fi^ — plane (dotted 
region). Such a region is elongated around the Q,o = 1 line identifying a flat geometry, Euclidean 
Universe. The blue lines define the age of the Universe for the considered models. The green shaded 
region is consistent (95% confidence contour) with the recent results of the high-redshift supernovae 

40 41 

surveys ' . 
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